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Let R be a simply connected region in En with smooth bounding surface S. For u suffi- 
ciently differentiable set of functions we derive a class of quadratic integral identities relating 
surface integrals of derivatives to integrals over R. These identities are a generalization of 
a first order identity given by L. Hormander (Compt. Rend. Douzieme Congr. des Mathe- 
maticiens Scandinaves Tenu a Lund, 1953, pp. 105-115) and L. E. Payne and II. F. Wein- 
berger (Pacific J. Math. (1958) pp. 551-573). As an example of an application of these 
identities we consider a solution u of the boundary value problem Au— pu=F in R and 
u=j on S. Here A denotes the Laplace operator and 0<p(x). We obtain pointwise a 
priori bounds for the derivatives of u in R in terms of a quadratic functional of an arbitrary 
function. Hence the Rayleigh-Ritz procedure can be used to make the error arbitrarily 
small. 



1. Introduction 

In a recent paper [6] 3 L. E. Payne and H. F. 
Weinberger give a method for obtaining bounds for 
solutions of second order elliptic boundary value 
problems. In that paper they give a generalization 
(see eq 2.4) of an integral identity of Rellich [7] which 
was essential to their method. This generalization, 
which was first obtained and applied to hyperbolic 
operators by L. Hormander [4], displays the highest 
derivatives in the form of a second order operator 
Lu^ia^u,)),^ This fact makes the identity useful 
in treating boundary problems for the corresponding 
differential equation Lu=F. 

Use was also made of this identity by Hubbard 
[5] to obtain bounds for membrane eigenvalues by 
finite difference methods. 

In section 2 we give a further generalization of 
Hormander's result, which involves higher deriva- 
tives on the boundary and where the highest order 
terms enter as derivatives of Lu. Section 3 gives 
an application of these higher order identities in 
obtaining pointwise bounds for the derivatives of a 
function u. The particular problem treated assumes 
a knowledge of Au—pu in a region R and u on the 
boundary C where A is Laplace operator and p^>0 
in jR+C. Because of the important physical appli- 
cation of thisjequation (see e.g. /Bergman and Schiffer 
[1]), this problem was chosen as an example to illus- 
trate a use of the higher order identity. More 
general equations could be treated with only tech- 
nical modifications. Other methods for obtaining 
explicit pointwise bounds for derivatives in such 
problems have been given by J. B. Diaz [3] and b}^ 
Payne and Weinberger [6]. In both cases the 
method given involves differentiation of the funda- 
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mental solution in Green's third identity, thus 
introducing more singular behavior of the bounds 
as the boundary is approached. In our paper the 
bounds for derivatives of all orders have the same 
behavior near the boundary, provided the boundary 
data is sufficiently differentiable. 

Extensions of the fundamental identity (2.4) can 
be derived for higher order operators by the same 
technique used here. One such, involving the 
Inharmonic operator A 2 , is applied in a forthcoming 
paper of J. II. Bramble and L E. Payne [2] to obtain 
explicit pointwise hounds in the Inharmonic Di- 
richlel problem. 

As has been previously pointed out in various 
places (see e.g. [(>]) bounds of the type obtained here 
can be used in conjunction with a Rayleigh-Ritz 
technique to approximate the unknown solution (or 
derivative) arbitrarily closely. 

2. A Class of Quadratic Integral Identities 



L 



L(u) = (a ii u, i ) {ii 



(2.1) 



be defined on a region R in E n with boundary C. 
The symbol ]f indicates covariant differentiation 
which coincides with fi (partial differentiation with 
respect to x 1 ) in Cartesian coordinates. A repeated 
index indicates summation from 1 to N. The sym- 
metric tensor a ij (x) is assumed to possess piecewise 
continuous derivatives of order M-\- 1 . In addition, 
let the eigenvalues of the matrix a ij be bounded 
away from zero and infinity in R. Hence there 
exists a positive constant, a, such that for all real 
nonzero (&, . . ., £„,) and all x in R f 



a' 1 2 &<*'%*,<* Z) &• 



(2.2) 



Under these conditions L is said to be uniformly 
elliptic. 
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Let f (x),A\ j (x), . . ., Al£(x) be symmetric tensors 
of the second order with piecewise continuous 
first derivatives in R. L. E. Payne and H. F. 
Weinberger have developed the following first order 
quadratic identity [6] and used it to obtain pointwise 
bounds in certain boundary problems. It is clear 
that 

(J m a i % iU { j ) [m = (f m a ij ) lm u {i u ]J + 2j m a l % t U\ mj 

= (f m a ij ) lm u ]i u u +2(f m a i % i u lm ) [j 

-2fT J a ij u [i u lm -2f m u lm L(u). (2.3) 

An application of the divergence theorem yields 

S{f m a i3 -2fa im }u {i u u v m ds=-2 \ f m u Vm L(u)dv 

JC J R 

+ f {{f m a t ') {m -2j\ m a im }u li u u dv. (2.4) 

J B 

In an analogous manner we can develop a homo- 
geneous quadratic integral identity involving second 
derivatives of u on the boundary C. Now 

(f m A^a k 'u lilc u ul ) im 

= (/M"a*0 ,„«, lk u ul +2f m A i 'a k % tk u Umt 
= +2(f m A<'a k \ ik u Un ) , l -2f m A l S i JLu) , , 
+ {(f'A i 'a kl ) lm -2(f'A i V m ) lm }u [/lc u u , 

+2f m A ii (aUu n ) ]( u Um . (2.5) 
After applying the divergence theorem we have 

&A tl {f*a t, -2f'a ,u } u, «tt|„iwfe 

-2 f /M < *tt, M (Lu), < efo+2 ( /'"^(</w u .) l? « l; „/fo 

J R J R 

+ f {(/M < V)i m -2(/'-A < ^* ra )i m }«i««i^e. (2.6) 

J R 

The corresponding identity which involves covariant 
derivatives of order M-\-l on the boundary is 

j> c AT . . . 4$*{/V'-2/a to >i« 1 . . . iM u, lh . . . }M v m dS 

= -2 f /M}* . . . A'^iLuh, . . . ««%, . . . ,„m* 
Jig 

+ f{(/ m ^r...^v') m 



-2Cf4 



* a *iii 



^M "<& ) I TO } 'Mj ii . . . itf frWl;, . . . injd® 



+ 2 f /^J* . . .Atf'iafMfou . . .^,* . . . JM mdv. 

J R 

(2.7) 

Since the integrands are in each case tensor in- 
variants, we may perform each integration in the 



most advantageous coordinate system. Since bound- 
ary conditions usually are given in terms of normal 
and tangential derivatives we shall display the deriv- 
atives appearing in the boundary integrals in such a 
form. Assume C to be a surface possessing M+l 
continuous derivatives as given in the parametric 
form 

x'=W, . . .,y N ~ l ), i=l, . . .,N. (2.8) 

The unit normal X*(y) satisfies the conditions 

G i J aX i =0 } o=l, . . ., JV-1, 

X i K i =\ (2.9) 

and the orientation of X i is taken inward. In what 
follows Greek indices will always range over 1, . . ., 
N—l, whereas Latin indices will refer to 1, . . ., N. 
We introduce geodesic normal coordinates (y l , . . ., 
y N ) in a strip immediately adjacent to the bounding 
surface. The transformation is given by 

x l =G\y\ . . ., yV-^+rXW, . . ., y N ~>). (2.10) 

The coordinate y N refers to a distance along the 
inward normal to C at the point (y . . ., y N ~ l ). 
The metric tensor has the form 

g«=(G£+y"X*)(G,l+y»X,l), 

9on= (G, k +y N X*)X*= hriX'X") ,«=0, 

<j NN =X k X k =-\ (2.11) 

The Christoffel symbols which involve the normal 
coordinate y N are 



\Nof 2 J by" 
\p«J ~2 



dy A 



{*}-{i/ w }-{^}-^ < 212 » 

The various second order covariant derivatives of u 
then have the form 



dy a dyP 



-{<£}**• 



1i",*jxr 



d 2 u 
llNa ~dy«dy N ~ 

d 2 u 



■1%}-' 



(dy»y 



(2.13) 



We indicate with a bar when the tensor is expressed 
in geodesic normal coordinates. Our operator takes 
the form 



Lu= (a iS u„) 1 1 = (5,*%,^ „. 



(2.14) 
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When L is the Laplace operator A we have 

Au=5 ij u, ij =u, NN +g a Pu la p, (2.1.5) 



since g aN =0 and g NN = 
takes the form 



1. Also tl 



lust 



Igrad u\ 2 =(u, N ) 2 +g a Pu, a u,p. 



invariant 



(2.16) 



We now express the boundary integral appearing in 
(2.4) in geodesic normal coordinates. 

-<£ {J N E ii -2j'a m }u, i M,jdS=-(£ Q N {{a^u, a u,^) 

-a m (u, N y]-2(J-u, a )(a iN u, i }dS. (2.17) 
If L = A then a ii =(j il and (2.17) becomes 

-^Jf N [(g^u, a u, B )-(u, N y]-2(f^u, a )u, N }dS. 

(2. IS) 

The surface integral in (2.0) can be treated in a, 
similar manner to obtain 

-(£ 3" Q N [a^ ka ^-a NN ^ m ) (u [lN )] 

-2(J"u lka )(a iN u {il )}dS. (2.19) 
In particular if j l ,A kl is chosen on Cso thai 

(2.20 

then For L = A (2.19) takes the simple form 

-(£jg a *u lN ji iN ?-(u lN „y}dS. (2.2i) 

In view of (2.15) we have thereby isolated an in- 
teresting combination of mixed normal and tangen- 
tial second derivatives on C in terms of second order 
tangential derivatives on C and certain integrals 
over R. 

In fact, a specific mixed derivative, say u lNa , can 
be isolated in this manner. Let 

/"=1, /«=0, 
A™=l,A ij =0 i,J9*(x (2.22) 

Then for L = A (2.19) becomes 

-<£ {grtu lffa u {ff p-(u m )*}dS. (2.23) 

As_we shall see in the next section, these choices of 
f\A kl will enable us to obtain new pointwise bounds 
for derivatives in certain boundary problems. 



The following notation will be used throughout 
this paper. Let/(x) be a piecewise continuous func- 
tion on R+C, then 



JM Z 



-m&x f(x) 

xeR+C ' 



J m =min f(x) 
xeR + C ' 



(2.24) 



For N=2j x i =O t (y 1 ) represents a plane closed 
curve. Equation (2.10) becomes 



x^&ffl+NnfiS) 



(2.25) 



where we have set y 1 =S (arclength along C) and 
y 2 =N (distance along the unit inward normal 
(n 1 , ri 2 ). If K{S) is the curvature of C then we have 



-(l-K(S)Ny 



g 1 



=<7 21 =0, g 22 =l. 



(2.26) 



The use of this coordinate system is, of course, 
restricted to a certain strip in the neighborhood of 
the boundary where in terms of it all points arc 
uniquely defined. 

3. Pointwise Bounds 

As was previously noted, Payne and Weinberger 
[()] made use of (2.4) in order to obtain bounds for 
solutions of second order elliptic partial differentia] 
equations. This identity enabled them to estimate 
the integral over a closed surface of the square of 
the normal derivative of a solution of a second order 
e(|uation in terms of integrals of the squares of the 
function and its tangential derivative over the 
surface. In order to bound the derivatives at. 
a point within the region they noted that one could 
differentiate the "Green's Identity, " being careful 
to define a "parametrix" in such a way that (his is 
permissible. The resulting estimates have the dis- 
advantage that the coefficients become infinite more 
rapidly as the boundary is approached than do 
those in the estimates for the value of the solution 
itself. This procedure has also been used by Diaz [3] 
in estimating derivatives at a point for solutions of 
Laplace's equation. 

One could proceed in an alternate manner. That 
is, instead of differentiating the "Green's Identity," 
simply write the identity using for the function to 
be evaluated at the point, the derivative of the 
solution. This leads to the estimation of surface 
integrals of squares of second derivatives of the 
function over the surface in terms of surface integrals 
of squares of tangential derivatives of the function 
over the surface. The identity (2.4) gives a means 
of obtaining the necessary inequalities. 

As an example of an application of these identities 
we consider the problem of obtaining a priori 
pointwise bounds for the derivatives of a function 
u for which Lu=Au—pu is known in It-\-C and u is 
known on C. Here II is a simply connected finite 
plane region bounded by the smooth closed curve 
C, A is the Laplace operator and p(x,y)>0 and 
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bounded together with its first derivatives in R+C. 
(The function p is taken to be positive in R+C for 
the sake of convenience. By a slight modification, 
similar results could be obtained for p>0.) We 
assume thai Lu and u are sufficiently smooth in 
R juk I on C respectively. 

As a starting point we use the "Green's Identity" 

«-J. r4 *H>.('s- r li)'* (31) 

where v is any sufficiently smooth function in B+(7, 
r is a point in the interior of R, T = — — lnr PQ (/> G = dis- 

tance from the point 7^ to another point (J), and — 

on 

is the outward normal derivative. Now let 
v=u,i=-^ — Then we have 

uAP^jjAu^lA+f^ g-r *"*)dS, (3.2) 

which may be written 

J R J R 



+* 






d//- 



W (3.3) 



We shall consider X,(P), since Xi(P) — u, t (P) is 
assumed known. Making use of Schwarz's inequality 
for vectors it follows that 

Xt(P)Xt(P) <K T Q jpuMA +D(u,u) 



du.f du.j 



IS 



where D(u,u) is Dirichlet integral, and 
K r =S f r 2 |giad A ^| 2 rf^+ f rydA 

J R J R 







(-3.4) 



Now we have 



k 



dn dn 



dS< 



(J) UtijUjijdS 



(3.5) 



(3.6) 



and in terms of the normal coordinate system intro- 
duced in section 2 (with the boundary coordinate 
taken to be arc length, i.e., y 1 = S ) y 2 =N), 

j)u, i/u, t /lS= (X)u] NN + 2u 2 NS +u 2 SJ3 dS. (3.7) 



Now on C, Av—pu=u lNN + u lss —pu, so that 
(j) u^jU^jdS^ (p [pu—Uiss— (Au—pu)] 2 

+u 2 ss dS+2(£ u 2 NS dS. 



Noting that 



U lS s= 



= ds 2 



K 



du d 2 u 



dN ds 



,+K 



du 
dn 



(3.8) 



(3.9) 



we have that 



+4lP M ^(^dS+2^uUsdS. (3.10) 

It follows from Green's Identity and the arithmetic- 
geometric mean inequality 



2aa- <a 2 a 2 -\ — p 



a>0 



(3.11) 



that 



if pu 2 dA+D(u,u)<6u^dS 



i r (Au—pu) 2 
+ 2j R ~y~ 



dA. (3.12) 



(V)mbining (3.10), (3.12) with (3.4) we obtain 



1 C (Au 
2J R 



i—puy 



dA 



Xi(p)Xi(p)<Kr 



+ 



(3.13) 

All terms on the right hand side of (3.13) excepl 
the last two are in terms of data. In order to bound 
the integral of the square of the normal derivative 
around the boundary we write the identity 

= ( \f, k S i3 -f j -f ti }u, 1 u,)dA-2l fu,Au<IA, 

J R J R 

(3.14) 

which has its left hand side essentially in the form 
of (2.18). Here a iJ =8 ij , n t and s t are the unit nor- 
mal and unit tangent vectors respectively and /* is 
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an arbitrary continuous, piecewise continuously dif- 
fereritiable vector field in ll-\~(\ We choose/* such 
that ( k n,c ^) on (\ (For example, if 11 is star-shaped 
with respect to the origin, then we may take /*=#*.) 
From (3.14) we obtain the inequality 



3 



(f 



■■■>sw 



dS 






]@y« 



+ «7 1 +X 1 +JBT a )D(tt,tt)- 



(rj'iov 



K { 



J r 



+^r f fPCAu-puydA (3.15 
&-2J r 



where Ci is any bound for the largest eigenvalue of 
the coefficient 'matrix [— f%8 ij +f,j+f,i], and .Ki and 
K 2 are positive. Let c=max [C^Xffl^M] and c hoose 

c+^+K^y. (Forexample, take 7^= ~ C+ ^ 2+4g 

and K 2 = ~~-) Then we have, using (3.12) and (3.11) 
with (3.15), 

| (/•»,),£ (g)' <*«<£[>», 

+ 4 ^f](£)'^+fJ/'«^-"»)^ 

+;j;, a ";"" s "- ii Si:,.";y'- < 3 - i6) 

or for simplicity 

+a 3 S (Au-pu¥dS (3.17) 



8c 2 






«3 



l/'. 



;} 



Finally, we need to obtain a bound for (Du 2 i N sdS. 

We make use of (2.23) and (2.6) assuming that the 
vector field j* and the tensor field A ij have been 
chosen as indicated. We then obtain the identity 

($(u 2 NS -u 2 ss )dS=-2 f j"A ij u ]jm {Au) u dA 

JC J R 

+ f i(f"A^) lm !f l -2(J l A^^u lik u ul dA. (3.18) 

J R 



We assume also that/' and A ij are bounded together 
with their derivatives in 11+0. (If the boundary 
is sufficiently smooth one can define these tensors 
(in the (S,N) system) as follows. Let 



f 0, ./' (\-KN) 
A 22 =A 12 =A 2l =0, A u (l-KN) 



for j\T<= 



for N>= 
K 



where if is a constant greater than the maximum 
curvature K M of C. 

Using (3.11) and Schwarz's inequahty we obtain 
from (3.18) 

CD u 2 NS dS<(p u 2 ss dS+ I (Au— pu),i(Au— pu),4. 1 

Jc JC J R 

+ | (pu)„(pu),M+ ( (" ik 'u Uk u ul dA (3.19) 

where 

C mi : -_2ffA mi A nj g mn \ (f\\<^) lm <f l -l(/A iJ ) lk . (3.20) 

The last term on the right hand side of (3.19) is 
bounded as follows: 

| C iJkl u {ik u ul dA<Bi u,, {J u, tj dA (3.21) 

where Bis a constant. After a moderate calculation, 

making use of the definitions of/' and A ij in the normal 
coordinate system, it is possible to obtain the bound 

B=l+K-K m +mnx(K/2\2(Kr/K T -K\ M ) 

where K' is the derivative of the curvature with 
respect to arc length. Now using the divergence 
theorem we have 

+ f (AufdA (3.22) 

J R 

where the n t are Cartesian components of the 
exterior unit normal. We may rewrite the boundary 
integral in terms of normal and tangential derivatives 
and obtain 

R U, i}U)i 4A-+2j cVnW% dS 

+ i K (iW + (W) **+ J>>^- (3 - 23) 
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Using (3.11) we have 

j R u, tJ u, tM <f c K (g) 2 dS+f c (gj dS 
+<j> (1+K) Q£f dS+2 £ (Au-jmydA 

+2 f {puYdA. (3.24) 

J R 

Using (3.11) and (3.12) we have the inequality 
f {pu)*dA<p M -((£ uHS 



dAf- (3.25) 



Combining (3.17), (3.24) and (3.25) we have 
f u^u^jdAKbduhlS+bd (j^JdS 

+ h(S(~^) 2 dS+b A f (Au-pu) 2 dA 

where 

b 1 =[(l+K) M +2p M ]a 1 +2p M , b 2 = [(l+K) M + 2p M )a 2 
+K M , h = l and h=2+[(l+K) M +2p M ]a,+2p M /p m 

(3.20) 

Now in (3.19) we are left to consider the terms 



(J) ufssdS and (pu),i(pu),i 
J c J r 



dA. 



In order to bound the latter we can write 
(pu) ,i(pu) yi dA 

J R 

= (p,iP, i u 2 +2pu, i p, i u+p 2 u, i u, i )dA 

J R 

<jJ^f^+p*)Qpu*+u, iU ,?}dA (3.27) 
where (3.11) has been used with a =-\l - 

M v 

Now (3.27) and (3.12) yield 
(pu), t (pu),idA 

J R 



(Au—pu) 2 



+j R i^-^i 2 d Ay (3. 



28) 



and hence, using (3.17) we have 
f (pu) , t (pu) ,4A<Ci& u 2 dS 

J R JC 

+c 2 S(^ydS+c 3 (£(Au-puydS (3.29) 



w here 



and 






Finally we have from (3.9), (3.11), and (3.17) 
& u] ss dS < 2K 2 M a, 6 u 2 dS+2K 2 M a 2 (j> (^~fdS 

+2(£ (0JdS+2K 2 M(h j {Au-pu)HA. (3.30) 

Combining (3.19), (3.21), (3.20), (3.29) and (3.30) 
we have 

+<7 4 (Au— pu) 2 dA-\-d b (Au— pu), t (Au— pu), t dA 

J R J R 

(3.31) 
where 

d 1 =2K 2 M a 1 +Bb 1 +c 1} d 2 =2K 2 M a 2 +Bb 2 +c 2 , 

d s =2+Bh, d i =2K 2 M a 3 +Bb,+c 3 and d 5 = l. 

Now (3.17) and (3.31) may be inserted into (3.13) 
to yield the desired bound, 

*MP) <K T i ft <§ c U*dS+p 2 f c (^J dS 

+ ft f c (^Sy dS+ ^ L ^ U ~^ HA 

+j8 4 f (Au-pu), i (Au-pu), i dA \ (3.32) 

where the ft's are explicitly determined constants. 

It should be pointed out that for the sake of 
simplicity, no effort was made to obtain the "best" 
such bound. Better bounds could be obtained, for 
example, by leaving certain known quantities under 
the sign of integration (compare equations (3.16) 
and (3.17)) rather than replacing them by a maxi- 
mum value. 
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4. Higher Derivatives in N Dimensions 

The process described in section 3 can be general- 
ized easily to higher derivatives. In bounding the 
Mth derivative we must require the existence and 
integrability of Mth derivatives of Au—pu, pu, and 
M+l tangential derivatives of u on C. We also 
assume that bounds have been obtained for the 
various lower order differential invariants. 

The mode of generalization will become apparent 
from the case M=2 where we assume that bounds 
have been obtained for 

J pu 2 dv, D(u,u), u^jU.ijdv, (J) ^— dS, (J) u yiJ u yij dS 

(4.1) 

as Mas done explicitly in two dimensions in the 
previous section. 

If co n denotes the surface area of the unit sphere 
in A 7 dimensions then 

r - (Ar-iU- N > 2 (4 2) 

is a fundamental solution of Au = 0. We can write 
Xu(P)=u, ti (P)- r(Au-pu), tJ dv= T(pu) fij dv 

J R J R 

Since T 2 dv does not exist for JV>4 we cannot apply 

J R 

Schwarz's inequality, as was done in (3.4), for dimen- 
sions higher than three. 

For N>4 we apply Schwarz's inequality as follows 

( f Tudv) 2 < f r N ~ x Y 2 dv f r~ {N ~ l) uHv 

\J R / JR JR 

where the singularity of r is chosen as the origin. In 
order to bound r~ {N ~ l) v?dv we use "Green's Iden- 

J R 

tity" 

- f [A(r- {N - 3) )u 2 -r- {N - 3) Au 2 ]dv 

J R 

=Boundary Integrals. 
Now using (3.11) we can easily obtain 

f r -(iv-i) u 2 ffo <__^_ f r-^- s) (Au-pu) 2 dv 

+ Boundary Integrals. 

This same basic technique may be used to handle all 
terms in (4.3) for N>4. 
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The boundary integrals occurring in the above in- 
equalities are the same as those to be dealt with 
below. 

We shall show how the desired bounds may be 
obtained for N<3 since in higher dimensions only 
technical details enter. 

Applying Schwarz's inequality, as was done in 
(3.4), we see that 

X (P)X ? -/P)</vY-( J (ptt) lW (pu)„/to 

+§ c U M lS+f c ^^<ls} (4.4) 

where K v is given by (3.5). 
We see that 

f0w) ><i (pii),,/fe<3(«^) (pu 2 

J R \ P /M J R 

+12(p, i p, t ) M D(u,u)-\-3p 2 M \ u,iju, tj dv, (4.5) 

J R 

but by assumption each integral on the right side of 

(4.5) has a known bound. Similarly we assume that 

a bound is known for (Du S ijfii- f tjdS. 

It now remains for us to obtain bounds for the last 
term of (4.4). We see that 

=CP g afi (f & U\ ayN u^ N dS- ] r 2(f) g afi U\ aNN u lBNN dS 
Jc Jc 

+^(u lNNN ) 2 dS. (4.6) 

From (2.15) we see that 

u ]<fNN = (Au—pu) ]<T + (pu)ia—g a ^u la ^. (4.7) 
The first two expressions on the right are known and 



U\aPy~ 



r d 2 u fk\-i 



dy y l_dy a dyP 



where 



du\ k _ 



W 



-{ay} T ^-{^Y a ^ (48) 

{i} *» (4 - 9) 



Um+- 



The first term on the right side of (4.8) is known. We 
can use the inequalities 



b iik %i^i k i<r 2 g ij g kl u liIC U\j h 



(4.10) 



where n, r 2 are upper bounds for the largest eigen- 
value in each case, to complete the bound for 



i 



g a ^u laNN u mN dS 



by a function of the quantities (4.1). 
Furthermore, from (2.15) we see that 

u WNN =(Au— pu) lN +(fm) lN — g^u ]am (4.11) 

and hence we can bound the last term of (4.6) in 
terms of known quantities and a bound for 



i 



g a Pg n u layN u { p SN dS. 



From (2.7) for M=2 and the choice 

A<#=g<*, A ia =A NN =0 

B™=1, B ij =0 ij?*N, (4.12) 

we have 

(V [g aB g y % ayN u [m —g^u {aNN u mN ]dS 

= - Jj m A ij B" l [(Au-pu)+pu] lik u ljlm dv 

+ f {(rA^B kl n [m -2(fA^B kl ) lr }u li!cr u lJls dV. 

J R 

(4.13) 
Hence 

(p g af, g n u layN : u im dS< (J) g all u laiNN u mN dS 

+ {-A^B kl [(Au-pu) lik (Au-pu) ul 

J n 

+ (pu) , ik (pu) | n ] dv+I(u { tjk ) , (4.14) 
where 

I(u Ujk )^jj(f m A i ^ l n ]m -2(fA i ^ l ) lr 

+ffA i *B kl }u {ikr u ljls dv. (4.15) 



But 



I(v> i tit) < ^s u ijJc u ij]c dv 

J R 

=tA Q) UijUjjicVkds— J (Au) tf jUijdv > 



4) g ii g kl u um u [jlN ds (4.16) 



+ e' 



where r 3 is again an upper bound for the largest 
eigenvalue of the coefficient matrix in (4.15). We 
can substitute (4.1) into (4.4) for e<l and achieve 
the desired bound for 



g^g lb u\ ayN u mN ds. 
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